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NOMENCLATURE 


e  =  measure  of  non-parallelism 

Independent  Coordinates 
t  =  time 

x  =  horizontal  spatial  coordinate 
y  =  vertical  spatial  coordinate 
=  ex 

Mean  Flow  Variables 

U(x^,  y)  =  U0(x-j ,  y)  +  eU^(x^,  y)  =  component  of  velocity 
V(x^ ,  y)  =  eV-j  (x-j,  y)  =  y  component  of  velocity 
P(xr  y)  =  P0(x1)  +  e  P](x1,  y)  =  pressure 
y)  =  Rq(x^  ,  y)  +  e  Rj(xj,  y)  =  density 
a(x^,  y)  =  speed  of  sound 
cp  =  specific  heat 
Perturbation  Variables 

ufxj,  y,  t)  =  (a0(xj,  y)  +  ea]  y))e10  =  x  component  of  velocity 

v(xp  y,  t)  =  (bQ(x1,  y)  +  eb^  (x^  y))e19  =  y  component  of  velocity 

p(xj,  y,  t)  =  (cQ(x1,  y)  +  ec-j  (xr  y))e  =  density 

P(x1 ,  y,  t)  =  ( dQ ( x ]  *  y)  +  ed1  (x-j,  y ) ) e 1 0  =  pressure 

s^xl»  y»  t)  =  (f0(Xl,  y)  +  ef1  (xlt  y))ei9  =  entropy 
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Perturbation  Variables  (Cont'd) 


0  =  k(x.j )  -  ojt 


u  =  frequency 
k(x-j )  =  wave  function 


Others 


=  matrix  operator 
=  adjoint  matrix  operator 


=  Hankel  function 


=  intermediate  functions  of  k 


vector  of  perturbation  components 


=  shear  layer  deflection 


u(x1),  v(x^)=  intermediate  functions  of  x-j 


=  adjoint  function 


=  perturbation  quantity 
=  perturbation  amplitude 
=  matrix  transpose 
=  variable  for  flow  outside  cavity 

=  variable  for  flow  inside  cavity 

=  cavity  length 
=  cavity  depth 
=  free  stream  Mach  number 


■ 
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I.  INTRODUCTION 

The  problem  of  a  compressible  flow  over  an  open  cavity  is  one  of  great 
practical  importance.  Examples  of  such  flows  are  those  over  aircraft  wheel 
wells  and  weapons  bays.  Bartel  and  McAvoy  (1981)  experimentally  measured 
sound  pressure  levels  as  high  as  170  dB  in  a  weapons  bay  environment.  Levels 
this  high  can  lead  to  structural  failure  and  will  cause  extreme  personal 
discomfort.  Thus,  methods  of  reducing  sound  pressure  levels  demand  immediate 
attention.  The  first  step  in  this  reduction  is  to  obtain  a  useful  analytical 
model  of  the  compressible  flow  over  an  open  cavity. 

Open  cavity  flows  have  been  investigated  analytically  and  experimentally. 
However,  many  analytical  models  have  severe  limitations  on  their  range  of 
application  or  do  not  match  experimental  data.  Some  of  the  best  analytical 
models  are,  in  fact,  semi-empi rical . 

Karamcheti  (1955)  performed  one  of  the  first  experimental  investigations 
on  a  simulated  weapons  bay.  He  obtained  discrete  frequency  radiation  for  both 
subsonic  and  supersonic  mean  flows.  This  result  has  been  verified  by  other 
investigators  including  Gibson  (1958),  Spee  (1966),  East  (1966),  and  Smith 
and  Shaw  (1974),  among  others. 

The  first  significant  analtyical  study  of  open  cavity  flows  was  that  of 
Plumblee, Gibson,  and  Lassiter  (1962).  They  identified  acoustic  resonances 
of  the  cavity  and  speculated  that  the  observed  discrete  frequencies  are 
identical  to  the  resonant  cavity  frequencies.  They  also  suggested  that  the 
driving  mechanism  of  cavity  oscillations  is  the  turbulent  shear  layer  spanning 
the  cavity.  Their  result  is  in  conflict  with  the  experimental  observations  of 
Karamcheti  (1955)  and  Rossiter  (1964)  whose  experiments  indicate  that  larger 
pressure  oscillations  are  obtained  for  laminar  boundary  layers. 

Based  upon  his  own  experimental  observations,  Rossiter  (1964)  proposed 
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a  vortex  shedding  model  to  explain  the  generation  of  cavity  tones.  He  was 
then  able  to  derive  a  semi -empirical  formula  for  the  tone  frequencies.  In 
addition  to  Rossiter,  Smith  and  Shaw  (1974)  and  Bartel  and  McAvoy  (1981) 
have  developed  semi -empirical  models  of  open  cavity  flows  based  upon  their 
experimental  observations.  These  models  do  reflect  the  data  obtained  by  the 
respective  investigators,  but  are  limited  in  application.  For  example,  the 
validity  of  Rossiter' s  model  appears  to  be  limited  to  the  Mach  number  range 
0.4  <  M  <  1.2. 

Analytical  efforts  at  modelling  open  cavity  flows  have  increased  within 
the  past  15  years.  Covert  (1970)  classified  cavities  as  shallow  or  deep 
depending  upon  their  length  to  depth  ratio.  He  argued  that  vorticity  tends 
to  excite  the  modes  in  the  direction  of  the  greatest  physical  dimension. 

Bilanin  and  Covert  (1973)  developed  an  analytical  model  of  a  shallow 
cavity  using  an  acoustic  monopole  to  model  the  trailing  edge  behavior.  Their 
model  is  based  upon  a  feedback  mechanism  first  proposed  by  Rossiter  (1964) 
where  the  pressure  field  of  the  monopole  drives  the  shear  layer.  The  cavity 
floor  influences  the  shear  layer  by  changing  the  amplification  rate  of  the  shear 
layer  motion  and  thus,  changing  the  excitation  frequency.  The  model  predicts 
excitation  frequencies  for  M  >  1.  Block  (1976)  extended  the  work  of  Bilanin 
and  Covert  to  include  effects  of  length  to  depth  ratio. 

The  appropriate  trailing  edge  behavior  has  created  some  controversy. 
Heller  and  Bliss  (1975)  analyzed  the  wave  motion  of  a  thin  shear  layer  over 
a  boundary.  They  superimposed  solutions  for  waves  moving  between  a  shear 
layer  and  solid  boundary  to  approximate  an  oscillating  cavity.  They  modelled 
the  trailing  edge  behavior  with  a  moving  piston  to  simulate  an  .entrainment 
process.  Heller  and  Bliss  argue  that  while  the  feedback  mechanism  model  of 
Bilanin  and  Covert  will  predict  excitation  frequencies  it  will  not  predict 
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whether  they  will  occur.  They  believe  that  an  entrainment  model  is  superior 
as  mass  addition  and  subtraction  must  occur  at  the  trailing  edge  as  the  shear 
layer  deflects.  Results  from  their  model  compare  favorably  with  results  from 
their  water  table  tests,  but  they  are  unable  to  predict  discrete  frequencies. 
Heller  and  Bliss  also  examined  methods  of  oscillation  suppression  using  vortex 
spoilers. 

Tam  and  Block  (1978),  in  a  study  directed  toward  open  wheel  well  noise, 
analyzed  the  open  cavity  problem  using  a  feedback  model  significantly  differ¬ 
ent  from  that  of  Bilanin  and  Covert.  However,  their  model  is  consistent  with 
the  water  table  observations  of  Heller  and  Bliss.  The  work  of  Tam  and  Block 
is  important  because  it  was,  until  this  study,  the  only  analytical  cavity 
study  to  include  shear  layer  thickness  effects.  They  used  a  mean  velocity 
profile  for  the  shear  layer  the  same  as  that  of  a  two-dimensional  free  turbu¬ 
lent  mixing  layer  near  he  trailing  edge  of  a  thin  flat  plate  and  assumed  the 
shear  layer  to  be  of  constant  momentum  th'ickness.  Their  results  compare  well 
with  the  experimental  data  of  Rossiter.  Tam  and  Block  also  show  that  shear 
layer  instabilities  could  be  the  origin  of  acoustic  energy  which  produces 
cavity  resonances. 

The  above  survey  of  literature  led  the  author  to  believe  that  a  unified 
analytical  theory  of  open  cavity  flows  is  lacking.  Some  theories  apply  only 
over  certain  ranges  of  significant  parameters  while  others  ignore  shear  layer 
thickness  effects  which  have  been  found  to  be  extremely  important.  In  view  of 
the  work  of  Tam  and  Block  the  author  also  decided  that  a  more  detailed 
analysis  of  the  shear  layer  including  thickness  effects  is  necessary.  Indeed, 
Bartel  and  McAvoy  recommend  a  thorough  mathematical  description.^  the  un¬ 
steady  pressure  distribution  in  the  shear  layer.  Thus,  the  author  proposed 
[Kelly  (1982)]  an  investigation  of  open  cavity  flows  including  shear  layer 


thickness  effects.  This  report  is  to  present  the  results  of  this  investi¬ 
gation  . 

The  author  decided  that  the  work  of  Tam  and  Block  was  a  significant  con¬ 
tribution  to  the  analysis  of  the  shear  layer  spanning  an  open  cavity.  Their 
results  are  valid  for  a  wide  range  of  Mach  numbers  and  length  to  depth  ratios, 
In  addition,  their  work  was  the  first  to  include  shear  layer  thickness 
effects.  They  analyzed  a  shear  layer  with  a  constant  momentum  thickness, 
but  noted  tli  ;  due  to  entrainment  the  thickness  of  the  shear  layer  increases 
in  the  downstream  direction.  They  also  used  their  shear  layer  instability 
analysis  to  predict  excitation  frequencies  from  an  equation  derived  assuming 
a  thin  vortex  sheet.  These  equations  are  summarized  in  the  next  section. 

From  his  prelimnary  research  the  author  realized  that  a  mathematical 
investigation  of  the  shear  layer  including  all  desired  effects  would  be  very 
involved  and  would  exceed  the  scope  of  the  project.  Thus,  it  was  decided  to 
focus  on  a  specific  aspect  of  the  shear  layer  based  upon  Tam's  analysis.  The 
author  recognized  that  there  are  two  major  drawbacks  to  the  work  of  Tam  and 
Block:  (1)  The  shear  layer  thickness  was  assumed  constant  across  the 
length  of  the  cavity  and,  (2)  The  eigenvalue  relation  they  use  should  be 
modified  to  include  thickness  effects.  Actually,  the  modification  of  their 
eigenvalue  relation  depends  upon  reanalyzing  the  flow  above  and  below  the 
shear  layer  assuming  a  shear  layer  of  finite  thickness.  Thus,  it  was  decided 
to  modify  the  work  of  Tam  and  Block  to  handle  a  shear  layer  whose  thickness 
increases  in  the  downstream  direction. 

Section  II  presents  a  statement  of  the  problem  under  consideration  as 
well  as  an  overview  of  the  work  of  Tam  and  Block. 

Section  III  presents  the  mathematical  analysis  of  a  compressible  inviscid 
shear  layer  whose  thickness  increases  in  the  downstream  direction.  The 


analysis  is  based  upon  non-parallel  stability  theory  using  the  method  of 
multiple  scales. 

Section  IV  presents  a  scheme  for  numerical  implementation  of  the 
analysis  in  Section  III. 

Section  V  discusses  the  method  and  presents  recommendations  for 


further  study. 


II.  Problem  Formulation 

The  objective  of  the  research  is  to  provide  a  mathematical  formulation 
of  the  shear  layer  over  an  open  cavity  including  thickness  effects.  A  side 
objective  is  to  evaluate  the  feasibility  of  using  a  mathematical  analysis  of 
the  shear  layer  in  a  finite  element  formulation  to  analyze  the  flow  in  all 
regions  around  the  cavity. 

A  literature  survey  indicated  that  the  only  study  including  shear  layer 
thickness  effects  was  made  by  Tam  and  Block  (1978).  They  first  analyzed  a 
thin  shear  layer  and  obtained  a  frequency-wave  number  relationship.  This  re¬ 
lationship  was  used  to  predict,  from  a  condition  derived  from  a  trailing  edge 
boundary  condition,  the  discrete  shear  layer  oscillation  frequencies.  They 
then  analyzed  a  shear  layer  of  finite  thickness  in  the  same  manner.  Their 
assumptions  were: 

1.  The  flow  is  inviscid. 

2.  No  mean  flow  exists  within  the  cavity. 

3.  The  shear  layer  is  of  constant  momentum  thickness. 

4.  The  mean  flow  in  the  shear  layer  is  parallel. 

5.  A  feedback  mechanism  occurs  due  to  the  deflection  of  the  shear 
layer  into  the  cavity.  Reflected  waves  excite  the  shear  layer. 

Tam  and  Block  were  able  to  predict,  with  favorable  comparison  to  experimental 
data,  discrete  oscillation  frequencies.  The  author  decided  that  a  first  step 
toward  a  thorough  mathematical  analysis  of  the  shear  layer  would  be  to  improve 
the  analysis  of  Tam  and  Block  by  eliminating  assumptions  (3)  and  (4).  Indeed 
the  shear  layer  thickness  increases  downstream  due  to  entrainment.  Thus  a 
shear  layer  whose  thickness  increases  downstream  was  analyzed  using  the  work 
of  Tam  and  Block  as  a  guide.  The  essential  features  of  their  work  is  presented 
below.  The  current  research  is  described  in  Section  III. 


The  flow  around  the  cavity  is  divided  into  three  regions:  (1)  The  flow 
inside  the  cavity;  (2)  The  flow  outside  the  cavity;  and,  (3)  The  shear  layer 
between  these  two  regions.  The  equations  governing  the  flow  inside  and  out¬ 


side  the  cavities  are 
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Where  a  +  subscript  refers  to  quantities  outside  the  cavity  and  a  -  subscript 
refers  to  quantities  inside  the  cavity. 

The  boundary  conditions  at  y  =  0  are 


<  it +  u- 1 11  =  v- 

(2.5) 

|a=  v_  +  A«(x-s)e'iut 

(2.6) 

p  =  p_  +  B6  (x-jje"'"1 

(2.7) 

Where  n  is  the  displacement  of  the  shear  layer,  A  and  B  are  constants  and  5 
is  the  point  where  the  shear  layer  is  excited.  Tam  and  Block  use  a  periodic 
line  source  at  the  trailing  edge  to  simulate  the  trailing  edge  interaction 
process.  They  calculated  the  pressure  field  inside  and  outside  of  the  cavity 
due  to  the  line  source  and  the  resulting  loading  on  the  shear  layer  due  to 
the  differences  in  pressure  and  velocity  above  and  below  the  shear  layer.  They 
are  then  able  to  obtain  a  relationship  between  the  shear  layer  deflection  at 
the  trailing  edge  and  the  source  strength.  The  pressure  attains  a  maximum 


value  when  the  shear  layer  deflection  is  most  negative.  Using  this  they 
derive  a  condition  which  defines  the  shear  layer  frequencies.  That  is 
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The  frequency-wave  number  analysis  is  determined  from  hydrodynamic  stability 
theory.  Equation  (2.8)  is  then  used  to  determine  the  discrete  frequencies. 

Simply  stated,  the  problem  at  hand  is  to  develop  a  technique  to  predict 
the  discrete  shear  layer  oscillation  frequencies  for  a  shear  layer  whose 
thickness  varies  with  x.  The  analysis  is  based  upon  hydrodynamic,  stability 
theory.  A  weakly  non-parallel  mean  flow  is  perturbed.  The  method  of  multiple 
scales  is  used  to  predict  an  eigenvalue  problem  which  governs  the  wave  number 


frequency  relationship.  The  wave  number  is  a  function  of  x,  k(x).  The 
eigenvalue  problem  is  solved  numerically  by  finite-difference  methods.  Thus 
a  numerical  relationship  between  k(L)  and  u  is  determined.  Equation  (2.8) 
is  used  to  determine  the  discrete  shear  layer  oscillation  frequencies.  The 
method  provides  a  correction  to  the  technique  of  Tam  and  Block  to  account  for 
varying  shear  layer  thickness  effects.  A  discussion  follows  the  analysis  to 
indicate  modifications  that  should  be  made  to  the  procedure. 

The  assumptions  made  in  the  analysis  are: 

(1)  The  flow  is  inviscid. 

(2)  The  flow  is  two-dimensional. 

(3)  No  mean  flow  exists  in  the  cavity. 

(4)  Equation  (2.8)  is  valid  for  this  problem. 

(5)  The  specific  heat  is  constant  across  the  thickness  of  the  shear 


III.  Shear  Layer  Modelling 


The  shear  layer  spanning  the  cavity  is  modelled  by  a  two  dimensional 
inviscid  compressible  layer  which  grows  in  the  downstream  direction.  The 
stability  characteristics  of  the  shear  layer  are  considered  when  it  is  ex 
cited  by  a  reflected  wave  of  frequency  u>. 

The  equations  governing  a  two  dimensional  compressible  inviscid  flow 


ie.  +  p  M.  +  u  Ip  +  p  iv  3p.  =  0 

at  M  ax  ax  M  ay  ay 


(3.1) 


x-Momentum 


/3u  .  au  .  „  3u\  ,  ap  n 
’  fe  ‘S  5^  a1 


(3.2) 


y-Mo men  turn 


°  (it  * u  fr  * v  If'  *  !y  * 0 


Entropy  Production 


(3.3) 


+  u  as  +  9s  = 
at  ax  ay 


Thermodynamics  Relation  for  an  Ideal  Gas 


dp  =  a^  dp  +  ~~  ds 
P 

The  mean  flow  is  assumed  to  be  a  steady  two-dimensional  flow. 
Thus,  the  mean  flow  quantities  (U,  V,  P,  R,  S)  satisfy 


(3.4) 


(3.5) 


w  \rm  \ tt 


U  p-  +  V  —  =  o 
a  x  ay 


9  Ra6 

dP  =  a  d  R  +  —  dS 
CP 


(3.9) 


(3.10) 


The  flow  variables  in  the  shear  layer  are  assumed  to  be  the  mean 


flow  component  plus  a  perturbation  quantity.  Thus 


u(x,  y,  t)  =  U(x,  y)  +  u(x,  y,  t) 


v(x,  y,  t)  =  V(x,  y)  +  v(x,  y,  t) 


p ( x »  y,  t)  =  R(x,  y)  +  p(x,  y,  t) 


p(x,  y,  t)  =  P(x,  y)  +  p(x,  y,  t) 


s(x,  y,  t)  =  S(x,  y)  +  s(x,  y,  t) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


Where  a  quantity  with  a  tilde  is  a  perturbation  quantity  and  is  assumed 
to  be  small  in  comparison  with  the  mean  flow  quantities.  Equations  (3.11) 
(3.15)  are  substituted  into  equations  (3.1  )-  (3.5).  Equations  (3.6)  - 
(3.10)  are  used  to  subtract  out  the  basic  state.  The  resulting  equations 
are  linearized  by  neglecting  all  terms  which  are  nonlinear  in  the  pertur¬ 
bation  quantities.  The  resulting  equations  are: 


~  +  R  —  +  p  —  +  U  —  +  u  —  +  R  —  +  p  — 
3t  3x  y  3X  3X  3X  3y  M  3X 


+  yi2.+  v— =0 

ay  3y 


(3.16) 


x-Mo men  turn 


R-|x-+Rll|^-+Ru|^-+p  uf^+Rvf^+Rvf^ 
at  ax  ax  ax  ay  ay 

+  ;  V  +  li  =  o 


(3.17) 


y-  Momentum 


R  |V  t  R  g  |1  t  R  -  ||  +  -  u  |X  t  R  V  §  t  R  0 

+  p  K  —  +  ^  =  0  (3.18) 

p  ay  ay 

Entropy  Generation 

3s  +  u  M.  +  u  ii  +  v  +  v  —  =  0  (3.19) 

at  ax  ax  ay  ay 


Thermodynamics  Relation 


P 


R  a2  <R  +  1^> 


(3.20) 
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Equations  (3.16)  -  (3.20)  form  a  set  of  five  equations  for  the  five 
unknown  perturbation  quantities.  Note  that  the  mean  flow  is  non-parallel. 

Thus  the  coefficients  in  the  differential  equations  (3.16)  -  (3.20)  are 
functions  of  both  spatial  variables,  x  and  y.  Standard  parallel  stability 
theory  would  suggest  assuming 

u  =  u(y)e1(kx  "  ult)  (3.21) 

where  k  is  the  wave  number,  to  is  the  excitation  frequency  and  u(y)  is  the 
amplitude  function.  However,  an  assumption  of  the  form  (3.21)  is  not  appro¬ 
priate  in  this  problem  as  the  coefficients  are  dependent  on  x. 

Instead,  an  approximation  to  the  solution  of  (3.16)  -  (3.20)  is  obtained 
by  following  the  work  of  Saric  and  Nayfeh  (1975).  The  non-parallel  effects  are 
assumed  to  be  weak  and  can  be  represented  by  a  small  parameter  e.  That  is, 
whenever  x  appears  in  the  mean  flow  quantities  it  appears  in  combination  with 
e  as  ex.  In  addition  the  non-parallel  component  of  velocity  V  must  be  of 
order  e  while  all  other  mean  flow  quantities  are  of  order  1.  Thus  it  is 
consistent  to  introduce  a  slow  scale  defined  by 

x  =  ex  (3.22) 

l 

Then  all  mean  flow  quantities  are  written  as  functions  of  x^  and  y.  Substitu¬ 
tion  of  (3.22)  into  (3.6)  -  (3.9)  yields 


From  equations  (3.23)  -  (3.26)  it  is  indeed  apparant  that 


U  =  0(1),  R  =  0(1),  V  =  0(e),  =  0(e),  S  =  0(1) 

oy 

Thus  it  is  also  convenient  to  expand  the  mean  flow  quantities  in  power 


series  in  e.  To  this  end 

U(x  ,y)  =  U  (x  ,y)  +  eU  (x  ,y)  +  ...  (3.27) 

1  0  1  11 

V(x  ,y)  =  eV  (x  ,y)  +  ...  (3.28) 

l  l  l 

R(x  ,y)  =  R  (x  ,y)  +  eR  (x  ,y)  +  ...  (3.29) 

1  0  1  0  1 

P(x  ,y)  =  P  (x  )  +  eP  (x  ,y)  +  ...  (3.30) 

1  0  1  11 

S(x  ,y)  =  S  (x  ,y)  +  eS  (x  ,y)  +  ...  (3.31) 

1  0  1  11 


The  method  of  multiple  scales  [Nayfeh  (1982)]  is  used  to  modify  the  assumed 
forms  of  the  perturbation  quantities  from  those  of  standard  stability  theory. 
To  this  end  let 


u(x  ,y,t)  =  a(x  ,y) 
l  l 

v(x,y,t)  =  b(x  ,y)  ei0(xi’t} 

1 0 ( X ! ,t) 

(x  ,y,t)  =  c(x  ,y) 
l  l 

p(x  ,y,t)  =  d ( x  ,y)  e1e^Xl,t^ 

l  l 

ie(xi,t) 


(3.32) 

(3.33) 

(3.34) 

(3.35) 


s(x  ,y,t)  =  f(x  ,y)  e 


(3.36) 


0  =  k(x  )  -  cut 
1 


(3.37) 


It  is  noted  that  both  the  amplitude  and  the  wave  number  are  assumed  to  be 
functions  of  the  slow  scale  x  .  It  is  also  noted  that  the  chain  rule  can 

l 

be  used  to  express  derivatives  with  respect  to  x  as 


3_  =  3 _L  +  3_  li 

3X  3X  dx  30  3X 

1 


f-  =  e  +  k(x  ) 

3X  3X  vi  30 

1 


(3.38) 


For  example 


t  fj-  [afx^yje16]  +  k(x()  |j  [afx^yje16] 


c«lr+  1ka]e 


It  is  also  convenient  to  expand  the  amplitude  functions  in  power  series  of 


a(x  ,y)  =  a  (x  ,y)  +  ea  (x  ,y)  +  ... 
1  0  1  11 


(3.39) 


b(x  »y)  =  b  (x  ,y)  +  eb  (x  ,y)  +  .. 
1  0  1  11 


(3.40) 


c (x  ,y)  =  c  (x  ,y)  +  ec  (x  ,y)  +  ... 
1  0  1  11 


(3.41) 


d ( x  ,y)  =  d  (x  ,y)  +  ed  (x  ,y)  +  ... 
1  0  1  11 


(3.42) 


f(x  ,y)  =  f  (x  ,y)  +  ef  (x  ,y)  +  ... 


(3.43) 
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Equations  (3.39)  -  (3.43)  are  substituted  into  equations  (3.31)  -  (3.36). 

The  resulting  equations  are  substituted  into  equations  (3.16)  -  (3.20)  making 
use  of  equation  (3.28).  Coefficients  of  like  powers  of  e  are  collected  in 
each  equation  and  set  equal  to  zero  independently.  The  results  are  summarized 
below 

Order  e°: 

Mass 


3r  3b 

i(kU  — oj )  c  +  b  — +  iR  ka  +  R  — ^  =  0 
o  o  o  oo  o  sy 


(3.44) 


x-Momentum 

3  U 

i ( kU  -u))R  a  +  R  b  t-2-  +  ikd  =  0 

0  0  0  0  0  dy  0 


(3.45) 


y-Momentum 

3d 

i(kU  -uj)R  b  +  =  0 

0  0  0  3y 

Entropy 


(3.46) 


i (kU  -u)f  +  b  =  0 

0  0  0  3y 


(3.47) 


Thermodynamic  Relation 

d  •  R  a2  (  ^  +  ) 

oo  R„  cp 


(3.48) 


Order  e1: 
Mass 


=  -  a 


>)c 

+  b 

3R 

— 9-  + 

ikR  a 

+  R 

3bi 

1 

l 

ay 

0  1 

0 

ay 

3Ro 

-  U 

3C0 

ike  U 

-  V 

3C 

0 

3X1 

0 

3xi  " 

0  1 

l 

ay 

3an  3b„  3V, 

- - R  ika  -  R  t— ^  -  c  - — 

io  i  ay  o  ay 


(3.49) 


X- Momentum 


Thermodynamic  relation 


Equations  (3.44)-(3.48)  represent  a  set  of  equations  to  solve  for  the 
zeroth  order  perturbation  amplitudes  ao,b0,c0,d0,  and  fo.  The  boundary  condi 
tions  in  the  y  direction  are  simply  that  the  flow  quantities  must  be  bounded 
as  y  -*  ±®  .  The  equations  can  be  written  in  the  following  form 
[L]  CU  =  (3.54) 

where 

Cl+J  -  Jo 
bo 
c0 
d0 
eo 

and  [L]  is  a  matrix  of  linear  operators  given  by 
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One  notes  that  the  operations  involved  in  [Lj  are  simply  multiplications 
by  known  functions  and  differentiations  with  respect  to  y.  Thus  it  is 
possible  to  consider  equation  (3.54)  at  a  given  spanwise  location  independent 
of  every  other  spanwise  location.  This  is  a  result  of  the  weakly  non-parallel 
approximation.  Arbitrary  functions  of  x-j  will  be  involved  in  the  solution  of 
equation  (3.54).  Indeed  it  is  convenient  to  write 

U0]  =  A(x^  [J  (xlty)]  (3.55) 

o 

where  A(x^)  is  at  this  level,  an  arbitrary  function  of  x.  However,  since 
differentiations  in  [L]  are  with  respect  to  y  only 

[L3U0]  =  CL]{  A[i0]l  =  A[L][i0] 

Thus  equation  (3.54)  yields 

CL]  Ci0  (x  v)]  -  0  (3.56) 

Equation  (3.56)  is  a  homogeneous  system  of  equations  to  solve  for  the 
components  of  [<t>0  (x^  ,y)].  The  matrix  [L]  is  a  matrix  of  linear  operators  that 
acts  on  the  components  of  a  five  dimensional  vector  whose  components  are 
functions  of  y  and  transforms  it  into  a  new  five  dimensional  vector.  Thus  the 
vector  spaced  is  defined  as  all  five  dimensional  vectors  whose  components  are 
complex  functions  of  a  real  variable  y  whose  range  is  from  -•*  to  +■*.  An 
inner  product  can  be  defined  on  this  vector  space  according  to 

(Cu],  [v])  =  f°° [u]T[V]dy 


— oo 


(3.57) 


23 


for  all  [u],  [v]  z'[).  The  bar  in  equation  (3.57)  denotes  a  complex  conjugate. 

For  a  given  excitation  frequency  w  it  is  desired  to  find  the  amplitudes 

of  each  variable  as  well  as  the  spatial  wave  function  k(x^).  Equation  (3.5) 

is  a  homogeneous  system  to  solve  for  these  amplitudes  at  each  spatial 

location  with  the  local  value  of  k(xj)  as  a  parameter.  However  this  system 
will  have  a  non-trivial  solution  only  for  certain  values  of  k(x^).  The  values 
of  k(xj)  for  which  this  non-trivial  solution  occurs  can  be  determined 
numerical ly . 

Equations  (3.49)-(3.55)  represent  a  set  of  equations  to  solve  for  the 
first  order  perturbation  amplitudes  aj,  ,  Cj ^  dj,  e^,  and  fj.  These 
equations  can  be  written  in  the  form 

[L][<j>,(xi,y)]  =  [g(x1,y)]  (3.58) 

where 


and 


U.  (*  ,  y)]  = 
1  i 


[g(x]  ,  y)]  - 


9]  (x  ,y) 
g2(x1  ,y) 

93(Yy) 

94!x,  Of) 

95(xry) 
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The  operator  [L]  in  equation  (3.58)  is  the  same  operator  that  appears  in 
equation  (3.54).  Equation  (3.58)  is  a  non-homogeneous  system,  whereas 
equation  (3.54)  is  a  homogeneous  system.  Furthermore  non-trivial  solutions  of 
equation  (3.54)  exist.  Thus  the  Fredholm  Alternative  implies  that  a  solution 
to  equation  (3.58)  exists  if  and  only  if  the  non-homogeneous  terms  satisfy  a 
solvability  condition.  In  particular  one  can  show  that  the  non-homogeneous 
terms  must  be  orthogonal  to  all  non-trivial  solutions  of  the  corresponding 
adjoint  problem. 

The  adjoint  of  an  operator  [L]  with  respect  to  a  given  inner  product  is 
the  operator  [L*]  that 

([L][u],[v])  =  (Cu],[L*]Cv])  (3.64) 

for  all  [u],[v]  e°°([L] ) .  Equation  (3.64)  applied  to  the  inner  product  given 
in  equation  (3.57)  becomes 


/”  (CL][u]I[v]dy  =  /”  [u]T( [L*][v])dy 


be  arbitrary  elements  of^lA  Then  using  the  definition  of  [L]  from  equation 
(3.55)  the  left  hand  side  of  (3.65)  becomes 


C  ([L][u])  [v]dy  =  OCIRqU,  +  ^  u2 


Integrating 

oo 

/ 


+  R0  3y  +  KW0  "  CKkU0  ”  w)  Rg^i 

3U0 

+  R0  3y"  U2  +  iku43  v2  +  C 1 ( ku0  •  “)  u2 
3u4  -  3Sn 

+  W  )]  V3  *  [  5TU2  *  i(kU0  -  u5  ]  v4 

2  V 

+  [-a  u3  +  u4  -  —  u5  ]  *5>  <ly 


by  parts  and  rearranging  terms  in  equation  (3.66)  yields 


(CL][u])T[v]dy  =  /*_  {[1R0;, 

3R 

+  1  ^ kuo  •  u)V2]ui  +  C  aT^i 

a  -  3Un 

•  fy  (R0V1}  'wV  i(kU0  '  w)  ^3 

+  aT  *4  ]u2  +  [i(kU0  •  “>  -  a%  3  u3 

9  V 

+  [1k;2  -  aT  +  ;5  ]  u4 


+  [i(kUn  -  uj)  v.  -  ~r —  vc  ]  uc }  dy 


i 


The  right  hand  side  of  equation  (3.67)  can  be  written  in  the  same  form  as 
the  right  hand  side  of  equation  (3.65)  with 


LL*]  - 


-i(kU  - 


9(Rq) 


■i(kUQ  -  u>) 


- i ( kU  -  o»)R  tt-2- 

'  o  '  o  3y 


-i  (kllg-O)) 


Let  [  <j>  *  ]  be  the  solution  to 
o 

CL*][i*]  =  0  (3. 

subject  to  the  requirement  that  all  components  of  [<{>*]  are  bounded  as  y  t 
The  Fredholm  Alternative  states  that  for  a  solution  of  equation  (3.58)  to 
exist  [g]  must  be  orthogonal  to  [$*].  That  is 


(Cg].[#)  =  o 


C  [greedy 


Equation  (3.70)  becomes,  upon  substitution  of  all  appropriate  quantities 


.  3Rn  3Cn 
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=  0 

(3.71) 


Equation  (3.71)  can  be  written  as 


P1  ^X1 )  dx-|  ‘  1y2^xl  ^  A  ^X1 )  =  0 


where  ^(x.j)  and  £(x^)  can  *5e  determined  by  comparing  equation 
equation  (3.71).  The  solution  of  equation  (3.72)  is 

Mp(x1  ) 

+i  !  -  7  dx.  .  .  ,  v 

U  (X,)  1  +lv(x,) 

A(x-| )  =  Aq  e  =  AQe 

where  Ag  is  an  arbitrary  constant  of  integration.  The  solution 
order  perturbation  quantities  can  be  written  as 


u(x] ,y)  =  AQe 


i(k(x1)  +  v(x1)  -  wt)a0(xry) 


i(k(x,)  +  v(x1) -ut)bnf  » 

v(xry)-A0e  1  1  0(xl’y) 


+  .  .  . 


p(x] ,y) 


p(xry) 


s(x] ,y) 


1(«x,)*v(x,)-.t)c0(Xi>y) 
+  .  .  . 

1(k(x,)tv(x1)..t)d0(X]>y) 
+  .  .  . 

f(k(.1)*.(«|)-.t)f0(  ) 

Ho 


(3.72) 
(3.72)  to 

(3.73) 

for  the  zeroth 

(3.74) 

(3.75) 

(3.76) 

(3.77) 

(3.78) 


The  sum 


X(X])  =  k(xn )  +  V  (X])  (3.79) 

represents  a  corrected  value  for  the  wave  number  to  include  the  effects 
of  the  shear  layer  thickness  increasing  downstream.  The  value  of  A ( L ) 
is  dependent  upon  the  frequency  at  which  the  shear  layer  is  excieted. 

For  a  given  w  the  procedure  described  in  the  next  section  can  be  used  to 
calculate  A (L ) .  A  numerically  defined  function  Q  (u>)  =  A (L )  for  a  given 
a)  is  defined.  This  function  can  then  be  used  in  equation  (2.8)  as  a 
correction  to  the  k(w)  used  by  Tam  and  Block. 


IV.  NUMERICAL  SOLUTION  OF  SHEAR  LAYER  EQUATION 


The  use  of  the  theory  described  in  Section  III  depends  upon  accurate 
numerical  evaluation  of  (a)  the  eigenvalue  k(xx ) ;  (b)  the  eigenfunction 

*  'ff 

[♦  (x  )];  (c)  the  adjoint  eigenfunction  [<i>0(x1)];  and,  (d)  the  wave  function 
correction  v(x1).  Each  of  these  quantities  must  be  evaluated  at  each  spatial 
location  along  the  length  of  the  cavity.  This  section  details  the  numerical 
solution  for  each  of  these  quantities. 

First  the  determination  of  the  eigenvalue  k(xJ)  and  the  eigenfunction 
[4>0(x1,y)]  is  considered.  These  quantities  are  obtained  by  solving  equation 
(3.54)  with  appropriate  boundary  conditions  at  each  spatial  location  such  that 
a  non-trivial  solution  is  obtained.  All  mean  flow  quantities  are  assumed  to 
be  known  before  the  solution  begins.  The  system,  when  written  in  matrix  form, 
leads  to  a  very  complicated  numerical  solution.  It  is  convenient  to  manipu¬ 
late  the  component  equations  of  equation  (3.54)  such  that  a  single  equation 
is  obtained  in  terms  of  one  variable.  The  component  equations  of  equation 
(3.54)  are 


i  (kU  -u>)cn  +  b 


3R  „  3b 

n  ^-+  iRnkan  +  R.  —  =  0 

0  3y  0  0 


0  3y 


(4.1) 


-  3U„ 

3d 
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i(ku0— )f.  *  b0^=  o 


2  .  C0  .  f0 


d0  =  ^0ai  (  R0  +  Cp  ^ 


Equation  (4.3)  is  solved  for  bQ  to  yield 

:  1  3do 

bo  =  (kU  -«)  3y~ 


(4.2) 

(4.3) 

(4.4) 

(4.5) 


(4.6) 
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Equation  (4.6)  is  substituted  into  equation  (4.2)  and  aQ  is  obtained  as 


ao  = 


R0(kU0-u»)'  [lkdo  +  kUp-o) 


9Un  3dn 

_ 0  _ o  1 

ay  ay  J 


(4.7) 


Equation  (4.6)  is  substituted  into  equation  (4.4)  to  yield 
1 


f  =  -  -  -  - 

0  _  /.  ..  .2 


3S„  3dn 
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Equation  (4.8)  is  substituted  into  equation  (4.5)  to  yield 


3  3d« 
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C°  =  7\(kU  -u,)2  ** 


(4.9) 


Equations  (4.6),  (4.7),  and  (4.9)  are  substituted  into  equation  (4.1)  to 
yield  the  following  differential  equation  for  dQ 
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(4.10) 


A  central  finite  difference  method  is  used  to  solve  equation  (4.10) 
for  the  eigenvalue  k(x:)  and  the  function  d^x^y)  at  discrete  values  of  Xj. 

The  range  in  the  y  direction  is  discretized  from  -HfySH  where  H  is  sufficiently 
large.  Convergence  tests  are  made  to  check  the  value  of  H.  Central  finite 
difference  formulas  are  used  to  approximate  the  derivatives  and  a  system  of 

A 

equations  is  derived  for  the  values  of  dQ  at  the  discrete  points.  The 
system  is  tridi agonal  thus  the  Thomas  algorithm  [Carnahan,  Luthur,  and  Wilkes 
(19)]  is  used.  Let  Ay  be  the  distance  between  points  in  the  y  direction  and 
define. 


(xt,  -  H  +  2Ay) 


(4.11) 


A  recurrence  relation  of  the  form 


dn .  =  G„  +  H„d  .. 
02,  2  2  02+1 


(4.12) 


is  assumed  where  additional  recurrence  relations  are  derived,  in  terms  of 


differential  equation  coefficients,  for  G  and  H  .  The  evaluation  of  the 

i  i 

G's  and  H's  occurs  in  a  forward  manner,  i  =  0,  ...,  N.  Then  the  d?'s  are 
evaluated  recursively  in  a  backward  fashion  i  =  N,  ...,  0  using  d^  =  0. 
However  the  G's  and  H's  are  in  terms  of  the  unknown  eigenvalue  k,  for  which 
an  initial  guess  must  be  made.  The  correct  solution  is  the  value  of  k  for 
which  the  differential  equation  and  its  boundary  conditions  has  a  non-trivial 
solution.  Values  of  G0  and  Hq  are  arbitrarily  selected  and  the  process 
begins.  The  value  of  dQQ  is  checked  against  zero.  If  it  is  not  close  enough 
to  zero  then  the  guess  for  k  must  be  changed.  A  Newton-Raphson  iteration  is 
used  to  iterate  the  values  of  k. 

The  remaining  perturbation  amplitudes  are  determined  using  equations 
(4.6)  -  (4.9)  and  numerical  differentiation. 

The  next  step  is  to  obtain  the  solution  of  the  adjoint  problem  given  the 
eigenvalue  k.  The  equations  defining  the  components  of  the  adjoint  vector  are 
obtained  from  equation  (3.68)  as 


-iRQa*  -  i(kUQ-io)R0b*  =  0 

(4.13) 

aRn  * 

a*  — o  _  (r  a*)  +  R  — -  b*  - 

3y  ay  v  o  '  o  ay 

as 

T(kU0-u,)c*  +  ^-d*  =  0 

(4.14) 

-i(kUQ-a))a*  -  a2f*  =  0 

(4.15) 

-i  kb*  -  ~  +  f*  =  0 

(4.16) 

R  a2 

-i(kUo-ui)d*  -  f*  =  0 

(4.17) 

Equation  (4.17)  is  used  to  solve  for  f*  in  terms  of  d*  as 
1CD(kU  -») 
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Equations  (4.15)  and  (4.18)  are  used  to  obtain 


a*  =  &s  d* 
Ko 


(4.19) 


Equations  (4.13)  and  (4.19)  are  used  to  obtain 


b*  = 


d* 


MkU<f“> 


(4.20) 


Equations  (4.16),  (4.18),  and  (4.20)  are  used  to  obtain 
iC, 
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(4.21) 


Equation  (4.14)  is  divided  by  (kU0-<u)  and  differentiated  with  respect  to  y. 
Equations  (4.19)  -  (4.21)  are  substituted  into  the  resulting  equation  to 
yield  the  following  differential  equation  for  d* 
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(4.22) 


Equation  (4.22)  can  be  solved  by  a  central  finite  difference  technique  once 
k  is  known.  Then  equations  (4.18)  -  (4.20)  are  used  to  solve  for  a*,  b*,  and 
f*.  Equation  (4.21)  along  with  a  numerical  integration  scheme  is  used  to 
solve  for  c*. 

The  process  outlined  above  is  continued  at  discrete  intervals  over  the 
entire  length  of  the  cavity. 

Numerical  integrations  in  the  y  direction  at  each  x^  location  are  per¬ 
formed  to  calculate  u^Xj)  and  y^Xj)  as  defined  from  equation  (3.71).  A 
numerical  integration  in  the  x ^  direction  is  then  performed  to  calculate  vfx^. 
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V.  DISCUSSION 


A  theory  has  been  presented  for  the  analysis  of  a  shear  layer  over  an 
open  cavity.  The  shear  layer  is  of  a  finite  thickness  which  increases 
downstream.  The  analysis  is  based  upon  non-parallel  stability  theory  using 
the  method  of  multiple  scales.  The  theory  is  used  to  predict  a  wave  number 
for  a  given  excitation  frequency.  The  wave  number-frequency  relationship  is 
then  used  to  predict  the  discrete  shear  layer  excitation  frequencies  using 
equations  developed  by  Tam  and  Block  (1978). 

A  numerical  scheme  for  implementing  the  theory  has  also  been  presented. 
Numerical  calculations  have  not  been  completed,  but  are  under  consideration, 
and  will  be  available  in  papers  published  from  this  work.  The  wave  numbers 
obtained  will  be  corrections  to  the  wave  numbers  obtained  by  Tam  and  Block  and 
will  be  compared  with  their  results. 

The  analysis  presented  is  not  intended  as  an  ultimate  analysis  of  the 
shear  layer  over  an  open  cavity.  The  author  recognizes  that  several  problems 
exist  with  the  analysis  and  many  questions  remain.  The  ultimate  goal  is  to 
attain  a  complete  mathematical  analysis  of  the  shear  layer  over  an  open 
cavity. 

One  question  concerning  the  analysis  presented  regards  the  mean  flow. 

Tam  and  Block  (1978)  note  that  the  mean  flow  in  the  shear  layer  has  never  been 
measured  such  that  the  results  are  useful  for  calculation  purposes.  Indeed 
they  use  a  standard  hyperbolic  tangent  velocity  profile  in  the  shear  layer. 

The  mean  flow  assumed  in  the  described  research  is  a  two  dimensional 
non-parallel  mean  flow.  The  non-parallelism  is  due  to  the  fact  that  the  shear 
layer  grows  downstream  due  to  entrainment.  However,  the  ratejjf  growth  of  the 
shear  layer  thickness  is  related  to  the  amplitude  of  the  cavity  oscillation 
and  is  not  predicted  by  the  analysis.  It  is  believed  that  a  thorough  analysis 


of  the  problem  will  yield  this  information.  Preliminary  ideas  are  presented 
below. 

The  stability  characteristics  of  the  shear  layer  are  only  required  over 
a  finite  distance  (0  5  x  <  L,  0  <  y  tL).  Thus  the  mean  flow  variables  can 
be  expanded  in  a  Taylor  Series  in  about  x^  =  0.  To  this  end 

8U0 

U0(xlty)  =  Uq(0 ,y )  +  x^  y^-(0,y)  +  .  .  .  (5.1) 

or 

U0(x1,y)  =  Aj  (y)  +  x1A2(y)  (5.2) 

The  nonlinear  terms  have  been  truncated  in  equation  (5.2).  The  other  mean 
flow  variables  can  be  written  as 


r0(xi  »y)  =  Bi  (y) 

+  x1B2(y) 

(5.3) 

Sg(x1  ,y)  -  c0(y) 

+  x-|C2(y) 

(5.4) 

V] (x1  ,y)  ■  D1 (y) 

+  x]02(y) 

(5.5) 

U] (x] ,y)  =  E1 (y) 

+  x]E2(y) 

(5.6) 

R] (*!  ,y)  =  F1 (y) 

+  x]F2(y) 

(5.7) 

S-|  (x1  ,y)  =  G1  (y) 

+  x1G2(y) 

(5.8) 

Equations  (5.2)  -  (5.8)  are  substituted  into  equations  (3.27)  -  (3.30)  and  the 
resulting  equations  substituted  into  equations  (3.23)  -  (3.26). 

A  simpler  version  of  the  mean  flow  analysis  is  to  assume  a  hyperbolic 


tangent  velocity  profile  whose  thickness  varies  with  x  .  To  this  end 


(5.9) 


U,  >  =  U  tanh(l  +  &  ■  \ 

(xlty)  -  e(x,J 

where 

0  (x-| )  =  Qq  +  x-|9 1  (5.10) 

The  result  of  the  method  described  is  the  shear  layer  stability 
characteristics  for  a  given  frequency  (i.e.,  a  frequency-wavenumber 
relationship.  This  relationship  is  then  used  to  predict  the  discrete  shear 
layer  frequencies.  The  author  feels  that  the  method  and  technique  are  good 
but  several  modifications  should  be  considered  before  the  method  is  applied 
for  design  analysis.  The  eigenvalue  prediction  equation  (2.8)  derived  by 
Tam  and  Block  is  based  upon  a  vortex  sheet  between  parallel  streams.  This 
prediction  equation  should  be  modified  to  include  shear  layer  thickness 
effects,  and  the  effect  of  shear  layer  growth.  This  will  involve  a  more 
detailed  mathematical  analysis  of  both  the  matching  between  the  shear  layer 
and  the  external  flow  and  the  fluid  behavior  at  the  trailing  edge  of  the 
cavity.  The  entire  fluid  mechanics  problem  must  be  considered  with  the  flow 
field  divided  into  three  regions.  The  non-parallelism  of  the  shear  layer  flow 
will  have  some  effects  on  the  external  flows  which  are  not  included  in  the 
eigenvalue  prediciton  equation.  Additionally,  this  equation  is  based  upon 
the  deflection  of  the  shear  layer  at  the  trailing  edge  of  the  cavity.  This 
concept  must  be  extended  to  handle  shear  layer  thickness  effects. 

The  boundary  condition  at  the  trailing  edge  must  be  considered.  The 
shear  layer,  while  being  excited  by  a  wave  of  frequency  -  may  either  reattach 
at  some  point  part  the  trailing  edge  of  the  cavity  or  deflect  into  the  cavity 
at  the  trailing  edge.  When  the  latter  occurs  due  to  entrainment  waves  are 
reflected  back  into  the  cavity  and  excite  the  shear  layer.  Thus  a  feedback 


process  is  born.  The  analysis  of  this  process  is  heavily  dependent  upon  the 
boundary  condition  at  the  trailing  edge.  An  analysis  of  the  shear  layer  with 
this  effect  should  be  performed.  Perhaps  it  will  be  discovered  that  a 
multiple  deck  like  analysis  is  necessary  where  higher  decks  are  needed  due  to 
the  trailing  edge  behavior.  Thus  a  complete  asymptotic  analysis  of  the  shear 
layer  should  be  performed. 

Some  investigators  have  attempted  to  solve  for  the  flow  variables  around 
an  open  cavity  using  finite  element  techniques.  It  has  been  suggested  [Wolfe 
(1982)]  that  a  mathematical  analysis  of  the  shear  layer  could  be  used  in  a 
finite  element  formulation  of  the  problem  to  predict  the  pressure  inside  the 
cavity.  This  analysis  is  the  first  step  toward  such  a  goal.  A  further  step 
in  this  process  is  to  conduct  an  analysis  of  the  shear  layer  with  this  goal  in 
mind.  That  is,  the  flow  variables  outside  the  shear  layer  will  be  solved  for 
by  a  finite  element  technique  while  the  variables  in  the  shear  layer  will  be 
provided  by  an  analytical  formulation.  However,  these  flow  regions  are 
intimately  connected  due  to  the  trailing  edge  interaction  process.  It  is 
possible  to  use  the  current  shear  layer  formulation  at  a  discrete  frequency. 
However,  the  amplitude  of  the  shear  layer  variables  is  not  currently  known  and 
must  be  solved  for  as  part  of  the  finite  element  routine. 

In  summary,  the  research  presented  is  just  a  first  step  toward 
understanding  the  shear  layer  over  an  open  cavity.  Much  remains  to  be  done, 
but  it  is  believed  that  future  work  is  manageable  and  will  lead  to  great 
benefits  to  reducing  sound  pressure  levels  in  cavities. 
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